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We search for the Zi(4050) + and ^2(4250) + states, reported by the Belle Collaboration, decaying 
to Xci7r + in the decays B° — > XciK~tt + and B + — > Xci-Ks 71 " 4 " where Xci — > J/4>~f- The data were 
collected with the BABAR detector at the SLAC PEP-II asymmetric-energy e + e~ collider operating 
at center-of-mass energy 10.58 GeV, and correspond to an integrated luminosity of 429 fb _1 . In this 
analysis, we model the background-subtracted, efficiency-corrected Xa 7r+ mass distribution using 
the Ktt mass distribution and the corresponding normalized Kir Legendre polynomial moments, 
and then test the need for the inclusion of resonant structures in the description of the Xci 7r+ 
mass distribution. No evidence is found for the -Zi(4050) + and Z2 (4250) + resonances, and 90% 
confidence level upper limits on the branching fractions are reported for the corresponding B-meson 
decay modes. 



PACS numbers: 12.39.Mk,12.40.Yx,13.25.Hw,14.40.-n 



I. INTRODUCTION 

The Belle Collaboration has reported the observa- 
tion of two resonance-like structures in the study of 
B° XciK-tt+ [1). These are labeled as Z x (4050)+' and 
Z 2 (4250)+, both decaying to xa^ + 0- The Belle Col- 
laboration also reported the observation of a resonance- 
like structure, Z(4430)+ -> vj(2S)tt + in the analysis of 
B —> ip(2S)Kir [H, These claims have generated a 
great deal of interest Such states must have a min- 
imum quark content ccdu, and thus would represent an 
unequivocal manifestation of four-quark meson states. 

The BaBAR Collaboration did not see the Z(4430)+ in 
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an analysis of the decay B —> vb(2S)Kir [6|. Points of 
discussion are: 

• The method of making slices of a three-body B 
decay Dalitz plot can produce peaks which may be 
due to interference effects, not resonances. 

• The angular structure of the B — > ijj(2S)Kir de- 
cay is rather complex and cannot be described ad- 
equately by only the two variables used in a simple 
Dalitz plot analysis. 

In the BaBAR analysis [f| , the B — > J/ip Ktt decay does 
not show evidence for resonances neither in J/ipn nor 
in J/ybK systems. All resonance activity seems confined 
to the Ktt system. It is also observed that the angu- 
lar distributions, expressed in terms of the Ktt Legendre 
polynomial moments, show strong similarities between 
B -> iP(2S)Ktt and B -> J/ipKn decays. Therefore, the 
angular information provided by the B — > J/i/jKtt decay 
can be used to describe the B —> iP(2S)Ktt decay. It 
is also observed that a localized structure in the ^(25)^ 
mass spectrum would yield high angular momentum Leg- 
endre polynomial moments in the Ktt system. Therefore, 
a good description of the yb(2S)TT data using only Kit mo- 
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merits up to L = 5 also suggests the absence of narrow 
resonant structure in the ip{2S)iT system. 

In this paper, we examine B — ¥ XciKtt decays follow- 
ing an analysis procedure similar to that used in Ref . @ . 
In contrast to the analysis of Ref we model the 
background-subtracted, efficiency-corrected Xci 7r+ mass 
distribution using the Ktt mass distribution and the 
corresponding normalized Ktt Legendre polynomial mo- 
ments, and then test the need for the inclusion of res- 
onant structures in the description of the Xci 7r+ mass 
distribution. 

This paper is organized as follows. A short description 
of the BABAR experiment is given in Sec. II, and the data 
selection is described in Sec. III. Section IV shows the 
data, while Sec. V and Sec. VI are devoted to the cal- 
culation of the efficiency and the extraction of branching 
fraction values, respectively. In Sec. VII we describe the 
fits to the Ktt mass spectra, and in Sec. VIII we show 
the Legendre polynomial moments. In Sec. IX we report 
the description of the Xci^ + mass spectra, while Sec. X 
is devoted to the calculation of limits on the production 
of the Zi(4050)+ and Z 2 (4250) + resonances. We sum- 
marize our results in Sec. XI. 



II. THE BABAR EXPERIMENT 

This analysis is based on a data sample of 429 ftV 1 
recorded at the T(45 l ) resonance by the BABAR detec- 
tor at the PEP-II asymmetric-energy e + e~ storage rings. 
The BABAR detector is described in detail elsewhere 
Charged particles are detected and their momenta mea- 
sured with a combination of a cylindrical drift chamber 
(DCH) and a silicon vertex tracker (SVT), both operat- 
ing within the 1.5 T magnetic field of a superconducting 
solenoid. Information from a ring-imaging Cherenkov de- 
tector is combined with specific ionization measurements 
from the SVT and DCH to identify charged kaon and 
pion candidates. Photon energy and position are mea- 
sured with a CsI(Tl) electromagnetic calorimeter (EMC), 
which is also used to identify electrons. The return yoke 
of the superconducting coil is instrumented with resistive 
plate chambers for the identification of muons. For the 
later part of the experiment the barrel-region chambers 
were replaced by limited streamer tubes [8(. 

III. DATA SELECTION 

We reconstruct events in the decay modes ||: 

B° -> XciK-tt+ , (1) 
B+ -> XciK° s 7T+ , (2) 

where Xci — ► J/07, and J/0 — ► P + H~ or J/0 —} e + e~ . 

For each candidate, we first reconstruct the J/ip by 
geometrically constraining an identified e + e _ or \x + \x~ 
pair of tracks to a common vertex point and requiring a 



X 2 fit probability greater than 0.1%. For J/ip — > e + e~ 
wc introduce bremsstrahlung energy-loss recovery. If an 
electron-associated photon cluster is found in the EMC, 
its three-momentum vector is incorporated into the cal- 
culation of m(e + e~) [lfj. The fit to the J/ip candidates 
includes the constraint to the nominal J/ip mass value 0] ■ 

A ifj candidate is formed by geometrically constrain- 
ing a pair of oppositely charged tracks to a common 
vertex (x 2 fit probability greater than 0.1%). For the 
two tracks the pion mass is assumed without particle- 
identification requirements. The fit includes the con- 
straint to the nominal mass value. 

The J/0, K , and 7r ± candidates forming a B meson 
decay candidate are geometrically constrained to a com- 
mon vertex and a % 2 fit probability greater than 0.1% is 
required. Particle identification is applied to both K and 
7r candidates. The K® flight length with respect to the 
B + vertex must be greater than 0.2 cm. 

A study of the scatter diagram E 1 vs. m(J/07) (not 
shown) reveals that no Xci signal is kinematically possible 
for E 1 < 190 MeV. Therefore, we consider only photons 
with a laboratory energy above this value. We select 
the Xci signal within ±2 a Xcl of the Xci mass, where 
o~ Xcl and the Xci mass are obtained from fits to the J/07 
mass spectra using a Gaussian function for the signal and 
a 2 nd -order polynomial for the background, separated by 
B and J/0 decay mode. The values of a Xcl range from 
14.6 MeV/c 2 to 17.6 MeV/c 2 . 

We further define B meson decay candidates using 
the energy difference AE = E* B — y/s/2 in the center- 
of-mass (cm.) frame and the beam-energy-substituted 
mass defined as m ES = yj ((s/2 + p 4 ■ p B )/Ei) 2 — p B , 
where (Ei,pi) is the initial state e + e~ four-momentum 
vector in the laboratory frame and y/s is the cm. en- 
ergy. In the above expressions E* B is the B meson candi- 
date energy in the cm. frame, and pb is its laboratory 
frame momentum. The B decay signal events are selected 
within ±2.0 (Jmss of the fitted central value, where the 
cmES va mes are listed in Table U and are determined by 
fits of a Gaussian function plus an ARGUS function [ll| 
to the data. 

The resulting AE distributions have been fitted with 
a linear background function and a signal Gaussian func- 
tion whose width values (ct^e) are a ^ so listed in Ta- 
ble HI Further background rejection is performed by se- 
lecting events within ±2.0 (r^E of zero. Table U also 
gives the values of event yield and purity, where the Pu- 
rity is defined as Signal/ '(Signal+ Background) . The AE 
distributions shown in Fig. [T] have been summed over the 
J/0 — s- and J/0 —> e + e~ decay modes. Clear 

signals of the B decay modes (1) and (2) can be seen. 
We obtain 1863 candidates for B° — > XciK~n + decays 
with 78% purity, and 628 B + — > XciK^tt + events with 
79% purity. A study of the AE and J/^7 spectra in the 
sideband regions does not show any B or Xci signal re- 
spectively. We conclude that the observed background is 
consistent with being entirely of combinatorial origin. 

The resulting J/07 invariant mass distributions for 
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FIG. 1: Distributions of AE for (a) B° -J- Xci^"7r + and (b) 
B + — > XciKsTT + summed over the J/tp decay modes; the Xci 
and tties selection criteria have been applied. The shaded 
areas indicate the signal regions. 



channels (1) and (2) are shown in Fig. [21 

In order to estimate the background contribution in 
the signal region, we define AE sideband regions in the 
intervals (7 — 9) <j/±e on both sides of zero. We obtain a 
"background-subtracted" distribution of events by sub- 
tracting the corresponding distribution for AE sideband 
events from that of events in the signal region. 



boundary. We observe a vertical band due to the pres- 
ence of the K* (892)° resonance and a weaker band due to 
the ^(l^SO) resonance. We do not observe significant 
accumulation of events in any horizontal band. 

The Dalitz plots for B + — > XciK®tt + candidates in the 
signal and sideband regions are shown in Fig. [4] and show 
features similar to those in Fig. [3] 



V. EFFICIENCY 

To compute the efficiency, signal MC events (full-MC) 
for the different channels have been generated using a 
detailed detector simulation where B mesons decay uni- 
formly in phase space. They are reconstructed and ana- 
lyzed in the same way as real events. We express the 
efficiency as a function of m(KTr) and cost9, the nor- 
malized dot-product between the \ci momentum and 
that of the kaon momentum, both in the Kir rest frame. 
To smooth statistical fluctuations, this efficiency is then 
parametrized as follows. 

We first fit the efficiency as a function of cos 9 in sepa- 
rate 50 MeV/c 2 intervals of m(Kw), in terms of Legendre 
polynomials up to L = 12: 
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FIG. 2: The J/tpj mass distribution for (a) B° -s> XciK~tt + 
and (b) B + — > Xci^s 71 "" 1 " candidates, summed over the J/ip 
decay modes. The m,ES and AE selection criteria have been 
applied. The shaded areas indicate the signal regions. 



For each value of L, we fit the ai(m) as a function of 
m{Kn) using a 6 th -order polynomial in m(Kir). The re- 
sulting fitted efficiency for B° decay is shown in Fig.(5ja). 
We observe a significant decrease in efficiency for cos 9 ~ 
+1 and 0.72 < m(if-7r+) < 0.92 GeV/c 2 , and for 

cos 9 1 and 0.97 < m{K--n + ) < 1.27 GeV/c 2 . The 

former is due to the failure to reconstruct pions with low 
momentum in the laboratory frame and the latter to a 
similar failure for kaons. A similar effect is observed in 
Fig. [SJb) for the the B + decay mode. 

In Fig. [5] we plot the efficiency projection as a function 
of m(xci^ + ) for channels (1) and (2), summed over the 
J/tp decay modes. We observe a loss in efficiency at the 
edges of the Xci 7r+ mass range. However these losses 
do not affect the regions of the reported Z resonances. 
Using these fitted functions we obtain efficiency-corrected 
distributions by weighting each event by the inverse of the 
efficiency at its (m(Kn), cos 9) location. 



IV. DALITZ PLOTS 

The Dalitz plots for B° — > Xci-^ _7r+ events in the 
signal and sideband regions are shown in Fig. [3J The 
shaded area defines the Dalitz plot boundary; it is ob- 
tained from a simple phase space Monte Carlo (MC) sim- 
ulation [l2| of B decays, smeared by the experimental 
resolution. For the sidebands, events can lie outside the 



VI. BRANCHING FRACTIONS 

We measure the branching fractions for B° —> 
XdK-n + and B+ -> XciK° s ir + relative to B° -> 
J/ipK + Tr~ and B + — > J / 'ip K^tt + , respectively. In this 
way several systematic uncertainties, (namely uncertain- 
ties on the number of BB mesons, particle identifica- 
tion, tracking efficiency, data-MC differences, secondary 
branching fractions) cancel. 
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TABLE I: Resolution parameter values from fits to the AE and toes distributions. 



Channel 






cr AK (MeV) 


am ES (MeV/c 2 ) 


events 


Purity % 


B°- 


> XclK 7T + 


(/< 


V) 


6.96 ± 0.34 


2.60 ± 0.10 


980 


79.3 ± 1.3 


B°- 




(e- 




7.81 ± 0.43 


2.77 ± 0.12 


883 


77.1 ± 1.4 


B+- 






6.65 ± 0.55 


2.65 ± 0.27 


299 


81.7 ± 2.2 


B+ - 




(e- 




7.52 ± 0.70 


2.65 ± 0.18 


329 


77.5 ± 2.3 
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FIG. 3: Dalitz plot for B° -> Xd^~ 

7r + in (a) the signal region and (b) the AE sidebands. The shaded area defines the Dalitz 

plot boundary. 



To obtain the yields, for each B decay mode we per- 
form new fits to the AE distributions using the full-MC 
lineshape for the signal and a linear background. The 
background-subtracted data are then integrated between 
±2.0 ^ ne correction for efficiency is obtained as 

described in Sec. V. A similar procedure is applied to 
the B° J/^K-n+ and B+ ->■ J/ipK°>rr + data. 

The branching fraction for Xci — > J/^l from Ref. @ is 
0.344 ± 0.015. Using this value, we obtain the following 
branching fraction ratios: 



B{B° -> XdK-ir+) 
B(B° -> J/tpK-n+) 



and 



B(B+ -> J/%PK°tt+) 



0.474 ±0.013 ±0.026, (4) 



0.501 ±0.024 ±0.028. (5) 



Systematic uncertainties are summarized in Table |TT] and 
have been evaluated as follows: 

1. We obtain the uncertainty on the background sub- 
traction by modifying the model used to fit the 
AE distributions. The signal was alternatively de- 
scribed by the sum of two Gaussian functions and 
the background was parametrized by a 2 nd -order 
polynomial. 

2. We compute the uncertainty on the efficiency 
by making use of the binned efficiency on the 



(m(2£V), cos 9) plane. In each cell we randomize 
the generated and reconstructed yields according 
to Poisson distributions. Deviations from the fitted 
efficiencies give the uncertainty on this quantity. 

3. We vary the bin size for the binned efficiency cal- 
culation. 

4. We include a systematic error due to the uncer- 
tainty on the Xci J 1^1 branching fraction Q. 

5. We assign a 1.8 % uncertainty to the 7 reconstruc- 
tion efficiency. 

6. We modify the AE and toes selection criteria and 
assign systematic uncertainties based on the varia- 
tion of the extracted branching fractions. 

We note that the systematic uncertainties are dominated 
by the uncertainty on the Xci J/ip"f branching frac- 
tion. 

The branching fractions measured in Ref. [|| are: 
B(B° -> J/ipK-TT + ) = (1.079 ± 0.011) x 10" 3 , (6) 



B(B+ -> J/i>K°ir + ) = (1.101 ± 0.021) x 10~ 3 , (7) 

where the latter value has been corrected for K® and 
K° ->■ 7r°7r° decays @. 
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FIG. 4: Dalitz plot for _B + — v XciKs^ in (a) the signal region and (b) the AE sidebands. The shaded area defines the Dalitz 
plot boundary. 
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FIG. 5: Fitted efficiency on the cos# vs. m(Kn) plane for (a) 

5° ->• Xci^"7T+ and (b) B+ -> X ciK° s -K + summed over the FIG - 6: Efficiency as a function of m( X cnv + ) for (a) B -> 
J/ip decay modes. XdK n + and (b) B + -> Xci-^s^ summed over the J/ip 

decay modes. 



TABLE II: Systematic uncertainties (%) for the B —¥ XciKir 
relative branching fraction measurements. 



Contribution B - 






1. Background subtraction 


1.6 


1.0 


2. Efficiency 


1.5 


1.6 


3. Efficiency binning 


1.1 


1.9 


4. Xci branching fraction 


4.4 


4.4 


5. 7 reconstruction 


1.8 


1.8 


6. AE and toes selections 


1.0 


1.0 


Total (%) 


5.4 


5.5 



B(B° -> Xd A'-7r+) = (3.83 ± 0.10 ± 0.39) x 10~ 4 . 

Multiplying the ratio in Eq. (5) by the B + — > 
J /ipK°ir + branching fraction in Eq. (7) we obtain 

B{B+ -> XciA°7r + ) = (5.52 ± 0.26 ± 0.31) x 10~ 4 , (9) 

so that, after all corrections, the branching fractions cor- 
responding to decay modes (1) and (2) are the same 
within uncertainties. 



Multiplying the ratio in Eq. (4) by the B a — s- 
J /ipK~ir + branching fraction in Eq. (6) we obtain 

B(B° -> xci K--K+) = (5.11 ± 0.14 ± 0.28) x 10~ 4 . (8) 

This may be compared to the Belle measurement [l|: 



VII. FITS TO THE Kir MASS SPECTRA 

We perform binned-^ 2 fits to the background- 
subtracted and efficiency-corrected Kn mass spectra in 
terms of S, P, and D wave amplitudes. The fitting func- 
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TABLE III: S, P, D wave fractions (in %), and x7 NDF 
(NDF = Number of Degrees of Freedom) from the fits to the 
Ktt mass spectra in B° — » XciK~ty + and B + — > Xci-K's 7r+ • 
The second P-wave entry in the two \ c \ channels corresponds 



to the fraction of A"* (1680). 


Channel S-wave 


P-wave 


P-wave 


X7NDF 


B° -> XciK-n+ 40.4 ± 2.2 
B + -¥ XciK° s % + 42.4 ± 3.5 


37.9 ± 1.3 

10.3 ± 1.5 
37.1 ± 3.2 

10.4 ± 2.5 


11.4 ± 2.0 
10.1 ± 3.1 


58/54 
55/54 



<j1000 

l 750 
O 

<N 500 



> 






1.2 

m(K~7T*) GeV/c 



m(K° s 7T + ) GeV/c 2 



tion is expressed as: 



r Gs(m) , p Gp(rn) , p Gp(m) 
J S J'Gs(m)dm T J P j G P (m)dm " t " ■'• D J'Gr>(m)dmJ ' 

where m = m(if7r), the integrals are over the full m(KTr) 
range, and the fractions / are such that 

fs + fp + f D = l. (11) 

The P- and D-wave intensities, Gp{m) and G.o(m), are 
expressed in terms of the squared moduli of relativis- 
tic Breit-Wigner functions with parameters fixed to the 
PDG values for if* (892) and if|(1430) respectively 0. 
For the S-wave contribution Gs (m) we make use of the 
LASS [01 parametrization described by Eqs. (11)-(16) 
of Ref.[| . 

The above model gives a good description of the data 
for the decays B — > JfyKn p|. However, for B — > XciKn 
the above resonances do not describe the high mass re- 
gion of the Kir mass spectra well. A better fit is obtained 
by including an additional incoherent spin-1 if * (1680) 0] 
resonance contribution. The fit results are shown by the 
solid curves in Fig. and the resulting intensity contri- 
butions are summarized in Table ILTTI In Figures [TJa) and 
[7]Jb) the contributions due to the if* (1680) amplitude are 
shown by the dashed curves. The XciKn decay modes 
differ from the corresponding J/ipKir and ijj(2S)Kir de- 
cay modes in that the .S-wave fraction is much larger in 
the former than in the latter. This was observed for the 
if* (892) region in a previous BaBAR analysis [Tij ]. 



VIII. THE Ktt LEGENDRE POLYNOMIAL 
MOMENTS 

We compute the efficiency-corrected Lcgcndre polyno- 
mial moments < Y® > in each Kit mass interval by cor- 
recting for efficiency, as explained in Sec. V, and then 
weighting each event by the Y®(cos6) functions. A sim- 
ilar procedure is performed for the AE sideband events, 
for which the distributions are subtracted from those in 
the signal region. We observe consistency between the 



FIG. 7: Fits to the background-subtracted and efficiency- 
corrected K-k mass spectra for (a) B° —¥ XciK~ty + and (b) 
P+ -> XoiK°n + . The K* (1680) contribution is shown in 
each figure by the dashed curve. 



B° and B + data. Therefore, in the following we combine 
the B° and B + distributions. 

This yields the background-subtracted and efficiency- 
corrected Legendre polynomial moments < Y£ >■ They 
are shown for L = 1, 6 in Fig. [8] We notice that the 

< Yq > moment is consistent with zero, as are higher 
moments (not shown). 

These moments can be ex pre ssed in terms of S-, P- 
and D-wave Ktt amplitudes [15(. The P- and D- waves 
can be present in three helicity states and, after inte- 
gration over the decay angles of the Xci, the relationship 
between the moments and the amplitudes is given by Eqs. 
(26)-(30) of Ref. [6j. We notice that, ignoring the pres- 
ence of resonances in the exotic charmonium channel, 
the equations involve seven amplitude magnitudes and 
six relative phase values, and so they cannot be solved in 
each m(K7r) interval. For this reason, it is not possible 
to extract the amplitude moduli and relative phase val- 
ues from Dalitz plot analyses of the if) Kir or XciKn final 
states. 

In Fig.[8]we observe the presence of the spin-1 if* (890) 
in the < Y® > moment and S-P interference in the 

< Yj > moment. We also observe evidence for the spin- 
2 if* (1430) resonance in the < F 4 ° > moment. There 
are some similarities between the moments of Fig. |8] and 
those from B — > J '/if) Ktt decays in Ref. @. However we 
also observe a significant structure around 1.7 GeV/c 2 in 

< Yj > which is absent in the B —> J/ipKir decays. We 
attribute this to the presence of the ^(1680) resonance 
produced in B — > XciKn but absent in B — > J/^pKir. 
The presence of scalar Z resonances should show up es- 
pecially in high < Y^ > moments. 

From the < Y® > we obtain the normalized moments 

< Y L N >= <^>, (12) 

where n is the number of events in the given m{Ki:) mass 
interval. 



11 



4\ M 



o 

X 



-0.5 





<Y 2 °> 


- \ 

It 




1 1 1 





O 0.5 1.2 1.9 0.5 

i 1 0.4 i— 



1.2 1.9 



> 
CD 



< 



Y 4 °> 



CO 

I -0.5 L 
0) 0.5 



"D 0.5 

CD 



^ h*r^Wr , . 



<Y 3 "> 



j'lf 1 Ttr w 



I I I I I I I _Q 2 I 1 

1.2 1.9 ' 0.5 

0.4 ,— 



1.2 



1.9 



O' 1 

"o 

5 



-0.5 





<Y 5 °> 






1 1 


I 1 



<Y B °> 




0.5 



_q 4 I 1 1 1 1 1 1 1 1 

1.9 ' 0.5 1.2 1.9 

„2 



rn(Kn) GeV/c 2 



FIG. 8: Legendre polynomial moments < Y° > for L = 1, 6 
as functions of Kir mass for B — > XciKn after background- 
subtraction and efficiency-correction. 



IX. MONTE CARLO SIMULATIONS 

We model B — > XciKn using the resonant structure 
obtained from the analysis of the K n mass spectra and 
Ktt Legendre polynomial moments. For this purpose we 
generate a large number of MC events according to the 
following procedure. 

• B — > XciKn events are generated uniformly in 
phase space [ll|. The B mass is generated as a 
Gaussian lineshape with parameters obtained from 
a fit to the data. 



TABLE IV: The value of x 2 / NDF for different MC-data com- 
parisons; "Y^ indicates the channel used to obtain the nor- 
malized moments. The "mixed" algorithm is explained in the 
text. The definition of window is given in Sec. X. 



Channel 








x7ndf 


B -t J/^K-k 


B - 


-> J/lpK 7T 


5 


162/152 


B -> XdKn 


B - 


-> XciKn 


5 


46/58 


B -> XciKn 


B - 


-> XciKty 


"mixed" 


63/58 


B — > XciKn window 


B - 


XciKiv 


5 


45/47 


B — > XclKty window 


B - 


-> XciKiv 


"mixed" 


56/47 



expression: 



w l= E <Y, N >Y l °(cos8). 

i=0 



(13) 



The moments correspond to the combined data 
from the decay modes of Eqs. (1) and (2). The 
< > are evaluated for the m(Kir) value by lin- 
ear interpolation between consecutive m{Kit) mass 
intervals. 



• The total weight is thus: 



(14) 



The generated distributions, weighted by the to- 
tal weight w, are then normalized to the number of 
data events obtained after background-subtraction and 
efficiency-correction. 

We first test the method using as control sample 
the combined data from B° — > J/tpK~n + and B + — » 
J/ip K !?7r + , where no resonant structure is observed in 
the J/-07T mass distributions jfjj]. In this case we gen- 
erate B — > J/i^K-k events and use the Kit resonant 
structure and Legendre polynomial information from the 
same channels. We compare the MC simulation to the 
J/ifiir + mass projection from data in Fig. We obtain 
X 2 /NDF = 223,162,180/152 for i max = 4,5,6 respec- 
tively. We conclude that L m ax = 5 gives the best de- 
scription of the data. 

We now perform a similar MC simulation for B —> 
XciKn using moments from the same channels. We ob- 
tain x 2 /NDF = 53,46,49/58 for L max = 4,5,6 respec- 
tively. The result of the simulation with L max = 5 is su- 
perimposed on the data in Fig.[T0l and the corresponding 
X 2 /NDF is given in Table ITVl The excellent description 
of the data indicates that the angular information from 



• We weight each event by a factor w m rjcjr) derived 
from the resonant structure in the Kit system de- 
scribed in Sec. VII (Eq. (10)), and displayed in 
Table HD 

• We incorporate the measured Kn angular structure 
by giving weight wj, to each event according to the 



the structures observed in the Xci n projection. This in- 
dicates the absence of significant structure in the exotic 
Xei7r + channel. 

We perform a MC simulation where, to the data from 
B° — > XciK~tt + , we add an arbitrary fraction (ss 25 %) of 
events which include a ^(4250)+ resonance decaying to 
XciTT- These Z 2 (4250) + events are obtained from phase- 
space MC B° — !• XciK ~7r + events weighted by a simple 
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FIG. 9: Background-subtracted and efficiency-corrected J/tpir 
mass distribution for the B — > J/ip Kit control sample with 
the superimposed curves resulting from the MC simulation 
described in the text. The solid curve is obtained using the 
total weight w obtained with L max = 5, the dotted curve by 
omitting the angular-dependence factor wl- 



FIG. 10: Background-subtracted and efficiency-corrected 
Xcitt mass distribution from B —¥ XciKtt. The solid curve 
results from the MC simulation described in the text, which 
uses of the moments from the same channels. The dotted 
curve shows the result of the simulation when the wl weight 
is removed. 



Breit-Wigner. We then compute Legendre polynomial 
moments for the total sample and use them to predict 
the Xci 71 " mass distribution as described above. The Xci^ 
mass spectrum for these events is shown in Fig. Illf a). 
We obtain X 2 /NDF = 103,91,88/58 for L max = 4,5,6 
respectively. Therefore, in the presence of a Z 2 (4250) + 
resonance, it is not possible to obtain a good description 
of the Xci 71 " mass distribution using L max = 5. We then 
increase the value of L max and obtain a good description 
of this MC simulation with L max = 15, as shown by the 
dashed curve in Fig. [Ilia) (x 2 /NDF = 57/58). 

We next test a "mixed" simulation where we use 
Lmax — 3 up to a Kn mass of 1.2 GeV/c 2 and L max = 4 
for the rest of the events. This choice is justified by 
the presence of spin and 1 resonances mostly in the 
low Ktt mass region, while the -ftT 2 (1430) contributes 
for m(Kir) > 1.2 GeV/c 2 . This simulation gives a 
satisfactory description of the B — > XciKir data with 
X 2 /NDF = 63/58 but gives a bad description of the MC 
sample of Fig. [□]>), yielding x 2 /NDF = 140/58. 

We now fit the MC sample including a simple Breit- 
Wigner (with the width fixed to the simulated value) to 
describe the Z 2 (4250)+ (Fig. Eft))). We obtain the solid 
curve, which has x 2 /NDF = 75/56, a Z 2 (4250)+ mass 
consistent with the generated value, and a yield consis- 



tent with the generated one. The dashed curve represents 
the background model from the "mixed" simulation. The 
MC test therefore validates the use of this background 
model for a quantitative evaluation of the upper limits 
described in Sec. X. 

The data-MC comparisons for the different simulations 
are summarized in Table IIVI 



X. SEARCH FOR Zi(4050)+ AND Z 2 (4250)+ 

We have shown, in the previous sections, that in the 
absence of Z resonances, the simulation with L max = 5 
gives a good description of the B —> J/ipKn and B — > 
XcI-^tt data. We now test the possible presence of the 
Zi(4050)+ and Z 2 (4250)+ resonances in B -> XciKir 
decay. Therefore we adopt the minimum L max config- 
uration ("mixed") described in Sec. IX and investigate 
whether something else is needed by the data. 

For this purpose we perform binned x 2 fits to the 
Xci7r + mass spectrum. In these fits the normalization 
of the background component is determined by the fit. 
We observe that this background model predicts an en- 
hancement in the mass region of the Z resonances. We 
then add, for the signal, relativistic spin-0 Breit-Wigner 
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TABLE V: Results of the fits to the Xciff mass spectra. N a 
and Fraction give, for each fit, the significance and the frac- 
tional contribution of the Z resonances. 



600 

ev 

O 
\ 

> 

4' 400 

O 
CM 

\ 

£ 200 

c 

o 
> 
o 



600 

in 
o 

\ 
> 

^ 400 

o 

CM 

\ 

i3 200 
c 
o 
> 

o 





4 4.5 

m(Xd tt) GeV/c 2 



FIG. 11: Background-subtracted and efficiency-corrected 
Xcitt mass distribution from B — > Xci-^Ti" which includes a 
simulated Z2(4250) + (vertical crosses). In (a) the distribu- 
tion with solid dots represents the B° —> XciK~iv + data 
component. The continuous curve is the result from the 
"mixed" simulation described in the text and obtained from 
the MC simulation. The dashed curve shows a simulation 
with L max = 15. (b) Result from the fit described in the 
text, which incorporates a Breit-Wigner lineshape describing 
the Z 2 (4250) + . The dashed curve represents the background 
model from the "mixed" simulation. 



functions with parameters fixed to the Belle values for 
the signals [l|. We compute statistical significance using 
the fitted fraction divided by its uncertainty. 

We first perform fits to the total mass spectrum. 

• Fit a) is shown in Fig. [T27a). and includes both 
Zi(4050)+ and Z 2 (4250)+ resonances. 



Data 



Resonance N a Fraction (%) x /NDF 



Total 



b) Total 



Zi(4050) + 
Z 2 (4250)+ 
Z(4150) + 



1.1 
2.0 
1.1 



1.6 ± 1.4 
4.8 ± 2.4 
4.0 ± 3.8 



57/57 
61/58 



c) Window Zi(4050)+ 1.2 3.5 ± 3.0 53/46 

Z 2 (4250)+ 1.3 6.7 ± 5.1 

d) Window £(4150)+ 1.7 13.7 ± 8.0 53/47 



• Fit b) is shown in Fig. IT27 b) , and includes a single 
broad Z(4150) + resonance. 

In both cases the fits give fractional contributions con- 
sistent with zero for the Z resonances. 

We next fit the Xci^ mass spectrum in the Dalitz 
plot region 1.0 < m 2 (K7r) < 1.75 GeV 2 /c 4 in order to 
make a direct comparison to the Belle results Fig- 
urcs[l2tc),(d) show the Xd 71 " mass spectrum for this mass 
region (labeled as "window" in Table |V| where the Belle 
data show the maximum of the reported resonance ac- 
tivity. This sample accounts for 25 % of our total data 
sample. Table ITVl gives the corresponding x 2 /NDF val- 
ues for the MC simulations described in Sec. IX, in this 
mass window. 

• Fit c) is shown in Fig. H"2T c). and includes both 
Zi(4050)+ and Z 2 (4250)+ resonances. 

• Fit d) is shown in Fig. [T2f d) , and includes a single 
broad Z(4150) + resonance. 

In each case the fit gives a Z resonance contribution con- 
sistent with zero. 

The results of the fits are summarized in Table [Vj and 
in every case the yield significance does not exceed 2a. 
Similar results are obtained when the resonance param- 
eters are varied within their statistical errors. 

We compute upper limits integrating the region of pos- 
itive branching fraction values for a Gaussian function 
having the above mean and a values, and obtain the fol- 
lowing 90% C.L. limits for the Zi(4050)+ and Z 2 (4250)+ 
resonances: 

B{B° -> Zi(4050)+if ") x £(Zi(4050)+ X d^ + ) (15) 



< 1.8 x 10" 
+ 



B(B" -> Z 2 (4250) + if ") x £(Z 2 (4250) + -> xa^ + ) (16) 

< 4.0 x 10~ 5 , 
B(B° -> Z+K-) x B(Z+ -> xci7r+) (17) 

< 4.7 x 10~ 5 . 

Systematic uncertainties related to the Z parameters 
have been ignored since they give negligible contribu- 
tions. The corresponding values for B + decay are «8% 
larger (see Eqs. (8) and (9)). 
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FIG. 12: (a),(b) Background-subtracted and efficiency-corrected XciT mass distribution for B — > Xa-K-K. (a) Fit with Zi(4050) + 
and Z2(4250) + resonances, (b) Fit with only the Z(4150) + resonance. (c),(d) Efficiency-corrected and background-subtracted 
Xoii" mass distribution for B — > in the Kn mass region 1.0 < m 2 (Kir) < 1.75 GeV 2 /c 4 . (c) Fit with Zi(4050) + and 

^2(4250) + resonances, (d) Fit with only the Z(4150) + resonance. In each fit the dashed curve shows the prediction from the 
"mixed" B — > \c\K-k simulation explained in the text. The dot-dashed curves indicate the fitted resonant contributions. 



Our measurements can be compared to the Belle re- 
sults 0: 



B(B° -» Zi(4050) + iO x £(Zi(4050)+ x c ivr+) (18) 

— ^ n +1 - 5 + 3 - 7 y m~ 5 

— °- u -0.8 -1.6 x 1U ' 

B(B° Z 2 (4250) + K-) x 6(Z 2 (4250)+ -S- x c itt + ) (19) 

= 4.0133 +w :l x io- 5 , 

(20) 



Given the large uncertainties, these branching fraction 
values are compatible with our upper-limit estimates. 



XI. CONCLUSIONS 

We use 429 fb _1 of data from the BAB4R experiment at 
SLAC to search for the Zi(4050)+ and Z 2 (4250)+ states 
decaying to Xci 7r+ m the decays B° — > Xci^ and 
B+ XdK°ir + , where Xci -> J/07 • 

We measure the following branching fractions for the 
decays B° — !• XciK~n + and B+ XciK°n + : 



B(B° -» X ciK-7r + ) = (5.11 ± 0.14 ± 0.28) x 10~ 4 , 
and 

B{B+ -> XciK°ir + ) = (5.52 ± 0.26 ± 0.31) x 10~ 4 . 
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In our search for the Z states, we first attempt to 
describe the data assuming that all resonant activity 
is concentrated in the Kit system. We use the decay 
B — > J/ipKTt as a control sample, since no resonant struc- 
ture has been observed in the J/ipn mass spectrum. In 
this case a good description of the data is obtained by 
a MC simulation which makes use of the known reso- 
nant structure in the Ktt mass spectrum together with a 
Legendre-polynomial description of the angular structure 
as a function of Kir mass. 

The same procedure is then applied to our data on 
the decays B — > xaKn and a good description of the 
Xci7r mass distribution is obtained. This indicates that 
no significant resonant structure is present in the Xci 71 " 
mass spectrum, as observed for the J/ip n mass distribu- 
tion [6|. We also observe that this background model 
predicts an enhancement in the mass region of the Z 
resonances. We then report 90% C.L. upper limits on 
possible B° — > Z + K~ decays. 

In conclusion, we find that it is possible to obtain a 
good description of our data without the need for addi- 
tional resonances in the Xci 7 ? system. 
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